Non-continuous Froude number scaling for the closure depth of a cylindrical cavity 
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A long, smooth cylinder is dragged through a water surface to create a cavity with an initially 
cylindrical shape. This surface void then collapses due to the hydrostatic pressure, leading to a rapid 
and axisymmetric pinch-off in a single point. Surprisingly, the depth at which this pinch-off takes 
place does not follow the expected Froude 1 / 3 power-law. Instead, it displays two distinct scaling 
regimes separated by discrete jumps, both in experiment and in numerical simulations (employing a 
boundary integral code). We quantitatively explain the above behavior as a capillary waves effect. 
These waves are created when the top of the cylinder passes the water surface. Our work thus gives 
further evidence for the non-universality of the void collapse. 

PACS numbers: 47.55.D-,47.55.db,47.11.Hj,47.35.Pq 



Many phenomena in fluid dynamics are known to be 
universal, allowing physicists to describe their final out- 
come without precise knowledge of the initial conditions. 
One prime example for such universality is the pinch-off 
of a liquid droplet in air [H, 0, 0] . The inverse problem, 
however, has recently turned out to be non-universal: 
When an air bubble pinches off inside a liquid, the local 
dynamics around the pinch-off point retains a memory 
of its creation until the very end [4|, [5|, [(| . In this Letter 
we show that also the global dynamics of this problem is 
characterized by the absence of universality. To this end 
we examine the air-filled cavity created when a solid ob- 
ject is rapidly submerged through a water surface. The 
walls of the cavity subsequently collapse due to hydro- 
static pressure from the liquid bulk. When the colliding 
walls meet, a violent jet shoots up into the air. Our ex- 
perimental and numerical evidence illustrates how uni- 
versality is broken through the interference of a second 
phenomenon unrelated to hydrostatic pressure. Surface 
waves created at the very beginning as the object passes 
the water surface exert a decisive influence until the end 
of the process. This leads to a discontinuous dependence 
of the pinch-off location on the velocity of the submerging 
object. 

In our experiment we drag a cylinder with radius 
Ro = 20 mm and length / = 147 mm through the surface 
of a large water tank using a linear motor connected to 
the cylinder bottom by a rod. With this setup we pre- 
scribe the cylinder velocity V between 0.5 and 2.5 m/s 
keeping it constant throughout the whole process. The 
relevant dimensionless parameter is the Froude number 
Fr = V 2 1 (Rod) with g = 9.81 m/s 2 , which in our experi- 
ment ranges between 1.2 and 32. 

The shape of the axisymmetric cavity is imaged with a 
high speed camera at up to 10,000 frames/sec, with the 
vertical coordinate z pointing upwards along the cylinder 
axis and r being the radial coordinate. Figure [T] shows a 
typical sequence of the cavity dynamics. We choose the 
starting position of the lower edge of the cylinder slightly 



below the water surface to suppress the splash (a). From 
there it is pulled downwards with high acceleration such 
that it has reached its prescribed speed before the top 
passes the water surface at t = 0. The submerging cylin- 
der creates an air-filled cavity whose side walls immedi- 
ately start to collapse (b). Magnification of the cavity 
walls as in (c) shows that they are not smooth surfaces, 
but exhibit pronounced ripples. The moment of pinch-off 
is shown in (d). 

From the high-speed video images we extract the clo- 
sure depth of the cavity z c for cylinder velocities up to 2.5 
m/s. For higher velocities the collapse becomes so violent 
that the video images do not allow for a clear determina- 
tion of the closure depth, thus limiting the range of our 
experiments. Plotting the closure depth over the Froude 
number as in Fig. [2j we find two asymptotic regimes. The 
regime for low Froude numbers obeys a scaling behavior 
z c ~ Fr a with an exponent a = ol\ « 0.1. This is in 
contrast to earlier experiments on the impact of similar 
cylindrical objects 0,11] • These showed a single continu- 
ous scaling behavior with a ~_Jr, for which a theoretical 
explanation could be given 

Only for high Froude 
numbers is our data consistent with this ^-scaling law. 
In between these two regimes we find a transition region 
without any experimentally discernable structure. 

To understand the underlying mechanism leading to 
this discontinuous behavior, we have conducted bound- 
ary integral simulations. This method [l(j assumes an 
inviscous and vorticity-free fluid. In fact, with the kine- 
matic viscosity v the global Reynolds number Re = 
RqV/v is of the order of 10 4 , while the local Reynolds 
number Re = RRjv defined with the cavity radius R for 
the point of minimum radius lies between 10 2 and 10 5 . 
We can thus safely neglect viscosity. Further, due to the 
small density of air as compared to water, its dynamic 
effects can be neglected for the study of the global cav- 
ity shape. Our numerical results can be compared to 
the high-speed videos without the use of any adjustable 
parameter. They are in good agreement with experimen- 
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FIG. 1: (color online) Snapshots of the cavity created by the submerging cylinder from the starting position (a) until pinch-off 
(d) taken at a cylinder speed of V = 1.5 m/s (Fr = 11.5). The area designated by the yellow rectangle in (b) is magnified in 
(c) illustrating the rippled surface of the cavity walls. Red and blue lines show the free surface and the cylinder, respectively, 
taken from the numerical simulation. In (a) and (b) video frames are connected to simulation time by matching the cylinder 
position. Due to the fact that the precise closure time sensitively depends on the amplitude of the capillary waves traveling 
over the cavity surface, for (d) it was more convenient to match the cavity closure times instead. 
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FIG. 2: (color online) The experimental closure depth as a 
function of the Froude number. The asymptotic regime for 
low Froude numbers scales with an exponent ai ~ 0.1 (red 
solid line). Only in the limit of high Froude numbers is the 
data consistent with an exponent of - (blue dashed line) as 
in @, 8, 9]. For the intermediate regime the experiments do 
not allow the extraction of a systematic behavior. 



tal observations as demonstrated by the colored lines in 

Fig. m 

These numerical simulations allow us to study the rip- 
ples observed in Fig.[T](c) in great detail. As the cylinder 
top passes the water surface, the rectangular corner be- 
tween the cylinder wall and the water surface is no longer 
held in place by the solid boundary of the cylinder. The 
newly created free surface thus possesses a corner with 
very high (initially infinite) curvature. Surface tension 
immediately tries to flatten this surface by pulling the 



corner diagonally inwards into the fluid bulk. This results 
in a shock similar to throwing a stone onto a lake, which 
consequently leads to the formation of capillary waves 
travelling out- and downward on the free surface. The 
downward waves can be observed in Fig. [1] (c). Such a 
shock creates a wave packet containing waves of different 
frequencies. Each of these waves spreads with a velocity 
c = ou /k given by the dispersion relation uj 2 — (cr/p) k 3 
of capillary waves, where uj is the angular frequency, k 
the wave vector, a the surface tension coefficient and p 
the density. Using the Fourier transform of this initially 
localized packet, we can apply stationary phase approx- 
imation and thus calculate the dominant wave vector /c* 
at a given distance x from the source at time t: 



Our simulations are detailed enough to allow for an 
accurate estimate of the dominant wave length A from 
the shape of the surface, cf. Fig. [3] (a). Comparison 
with Eq. (pp) as in the inset of Fig. [3] (a) confirms that 
the observed ripples on the surface are indeed capillary 
waves originating from the corner point as the cylinder 
top passes the water surface. 

The closing of the cavity is driven by hydrostatic pres- 
sure which acts on every point of the vertical free sur- 
face accelerating it inward as soon as the cylinder has 
passed. This accelerating force increases with the depth. 
Thus, points near the top surface start moving early with 
a small acceleration, while deeper points start with in- 
creasing delay, but higher acceleration 0, 0] . On a rip- 
pled surface this process favors the wave crests over the 
other points. The resulting closure depth will thus be 
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FIG. 3: (color online) (a) Capillary waves at t = 5.35 ms after the cylinder top has passed the water surface with V = 3 m/s 
(Fr = 45.9) as obtained from the simulation. Wave crests are marked by red circles and allow an estimate of the dominant 
wave length A = 27v/k at a given position. The inset compares these wave lengths for the downward (blue circles) and outward 
(black crosses) waves to the theoretically expected behavior for capillary waves of Eq. (pQ) (red line), (b) The natural cavity 
dynamics (red) as a superposition of two hypothetical settings: without surface tension (black) and without gravity (blue) for 
V = l m/s (Fr = 5.1). 



determined by a combination of (i) hydrostatic pressure 
induced by gravity and (ii) capillary waves created by 
surface tension. 

The numerically obtained closure depth shown in Fig.[H 
convincingly reproduces the experimental results, while 
for very high Froude numbers only numerical data is 
available. The first regime shows a scaling behavior with 
aq « 0.08, sufficiently close to the experimental value. 
The second regime scales with a 2 ~ 0-43, which is some- 
what off the expected value of | from @, H, @] , but well 
compatible with our experimental observations. The sim- 
ulations even make the identification of an intermediate 
regime possible, which due to its small range in Froude 
numbers cannot be clearly observed in experiments. The 
slight shift between numerical and experimental data on 
the Froude axis can be attributed to small contamina- 
tions which lower the surface tension of the water. 

The insets in Fig. [H elucidate precisely how the cap- 
illary waves lead to the discontinuous jumps between 
the different regimes of the closure depth: For Froude 
numbers near the regime transition three different local 
minima of the radius come very close to meeting their 
counterpart on the opposite side. In the first regime the 
uppermost of the three minima closes first and thus de- 
termines the closure depth, while in the second regime 
the one located in the middle is the fastest to reach the 
central axis. Finally, when the lowest minimum closes 
before the other two the third regime is attained. 

Since the behavior of capillary waves is determined by 
surface tension, one expects that modifications of the sur- 
face tension coefficient a should significantly alter the 
closure depth in the first and any intermediate regimes. 
For the last regime capillary waves are irrelevant and the 
scaling behavior of the closure depth can be derived inde- 
pendent of surface tension @, H, 0] . In numerical simula- 
tions a significant change of the surface tension coefficient 
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FIG. 4: (color online) Comparison of the experimental clo- 
sure depth (black crosses) with the numerical data (blue dia- 
monds). The dependence of the closure depth on the Froude 
number exhibits three clearly distinct regimes. The insets 
illustrate the shape of the cavity at pinch-off for one repre- 
sentative of each regime (axis are the same for all insets). The 
regimes are determined by which local minimum first meets 
its counterpart on the central axis. 



is feasible and we have performed simulations with a ten- 
fold increase and decrease of a as compared to the natural 
value of 72.8mN/m for water as well as for a hypotheti- 
cal liquid without any surface tension (a = 0). As Fig. [5] 
demonstrates, the Froude number ranges for the three 
regimes are found to depend indeed strongly on the value 
of the surface tension coefficient. The length of the first 
regime significantly enlargens for a higher surface tension 
coefficient. As expected, the last regime is almost unin- 
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FIG. 5: (color online) The numerical closure depth as a func- 
tion of the Froude number for different surface tension coef- 
ficients: (Ji = 728 mN/m (black circles), 02 = 72.8 mN/m 
(blue diamonds), (73 = 7.28 mN/m (red crosses) and a a — 
(green squares) . A higher a leads to a significant extension of 
the first regime. The inset shows the onset of the high Froude 
number regime as a function of the Bond number with the 
red line depicting the expected scaling law Fr tr ans ~ Bo _3//4 . 

fluenced by changes in surface tension. For low surface 
tension merely the two asymptotic regimes exist and the 
intermediate regime is not observed anymore. The lim- 
iting case completely without surface tension possesses 
only one regime, consistent with the Fr 1 ^ 3 -scaling. 

Using Eq. (pQ) with x = z c and k* ~ l/z c we can es- 
timate the relevant time scale for the capillary waves as 
t w ~ yj p/az^ 2 while the relevant time scale for cav- 
ity closure is t c ~ z c /V. The onset of the high Froude 
number regime is now readily found by equating the ra- 
tio t w /t c to a constant of order 1. Introducing the Bond 
number Bo = gR^p/cr and making use of the theoretically 
expected scaling z c ~ Fr 1 ^ 3 this yields Fr tra ns ~ Bo _3//4 
which is in good agreement with our observations as 
shown in the inset of Fig. [5l 

The cavity dynamics completely without surface ten- 
sion and thus deprived of all capillary waves is illustrated 
by the black lines in Fig. [3] (b). The blue lines depict the 
evolution of a free surface starting with a (nearly) rect- 
angular corner without gravity allowing us to study the 
formation of capillary waves in an isolated setting. The 
cavity dynamics under realistic conditions (g = 9.81 m/s 2 
and o~ = 72.8 mN/m) is shown in red and can clearly 
be identified as a superposition of the above mentioned 
limiting cases. In the first instants, the real dynamics 
is almost identical to the one without gravity with its 
main feature being the capillary waves. Later in the pro- 
cess hydrostatic pressure becomes more important until 
finally the cavity approaches the shape of the pure grav- 
ity simulation with the capillary waves superposed on the 



cavity walls. 

Finally, to examine the local dynamics around the 
pinch-off point we extract from the simulation the neck 
radius and the radius of curvature at closure depth 0, EH . 
Both scale with the time r = t c — t until pinch-off with 
scaling exponents a ne ck = 0.54 ... 0.57 and ^curvature = 
0.28 .. . 0.43, which agrees well with the observations for 
an impacting disk in [4]. The diffferent scaling behavior 
for the two length scales suggests that the local dynamics 
exhibits similar non-universal features as in [4|. 

In conclusion, we have shown that capillary waves cre- 
ated when a submerging object passes the water surface 
have a strong and lasting influence on the global dynam- 
ics of the cavity. This influence remains observable until 
the very end of the cavity collapse manifesting itself in 
clearly distinct regimes of the closure depth as a function 
of the submerging velocity. We have thus presented ev- 
idence that the pinch-off of air in water - in contrast to 
the pinch-off of a water droplet in air - is non-universal 
not only in a local 0,0, but also in a global sense. Since 
capillary waves are an unavoidable consequence of dis- 
turbances on a water surface, we expect that the effects 
elaborated in this Letter will be of relevance to a wide 
range of related phenomena. 
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